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128,130 Te and 15t) Nd nuclei for the + — ► + transition is studied in the PHFB model in conjunction with 
the pairing plus quadrupole-quadrupole plus hexadecapole-hexadecapole effective two-body interaction 
and the effect of the latter is investigated on the calculation of nuclear transition matrix elements M2v 
The reliability of the intrinsic wave functions of parent and daughter nuclei involved in the (/3~f3~ J 



decay of above mentioned nuclei is established by obtaining an overall agreement between a number of 
theoretically calculated spectroscopic properties, namely the yrast spectra, reduced B(E2:0 + — > 2+) 
transition probabilities, static quadrupole moments Q(2 + ) and g-factors 3(2+) and the available 
experimental data. The effect of deformation on Miv is also investigated to inveterate its inverse 
relation with nuclear deformation. 
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1. Introduction 

The nuclear fjf3 decay is a rare second order semi-leptonic transition between two even Z-even N 
isobars i|X and Z ± 2 Y involving strangeness conserving charged weak currents. The (3(3 decay can be 
broadly classified into four experimentally distinguishable modes, namely two neutrino double beta {(3(3) 2v 
decay p], neutrinoless double beta ((3(3) 0l , decay [2], single Majoron accompanied neutrinoless double 
beta {0(34>)q v decay [3] and double Majoron accompanied neutrinoless double beta {(3(3(f)(f)) 0l/ decay [4]. 
The(P(3) 2l/ decay conserves the lepton number L exactly and is an allowed process within the standard 
model of electroweak unification (SM). In ((3(3) 0u decay, the conservation of lepton number is violated by 
two units and it is possible in models beyond the SM, namely GUTs (left-right symmetric SO(10), E(6) 
etc.), Rp— conserving as well as violating SUSY models, leptoquark, compositeness and sterile neutrino 
scenarios. In fact, the experimental observation of ((3(3) Qu decay would immediately imply that neutrinos 
are Majorana particles and all the present experimental activities are directed towards the observation 
of this particular decay mode. The {(3(3(f) 0v and ((3(3<jxj)) 0ll decay modes are possible processes in nine 
Majoron models as discussed by Bamert and co-workers [5]. 

These decay modes can proceed via emission of two electrons ((3~ (3~), emission of two positrons 
((3 + (3 + ), electron-positron conversion (e/3 + ) and double electron capture (ee). The latter three are 
energetically competing modes and we refer to them as e + (3(3 modes. There are 35 (3~(3~ and 34 e + (3(3 
emitters. Further, the transitions in (3(3 decay modes may be from + — » J + states. Presently, the 
study of (3~ (3~ decay for the + — » + transition is the most preferable as the other decay rates are 
suppressed due to kinematic reasons. However, the study of other decay modes will be of importance 
in distinguishing the role of different mechanisms involved in (P(3) 0v decay once it is observed. Thus, 
the experimental as well as theoretical study of nuclear (3(3 decay is quite wide in scope and has been 
excellently reviewed over the past years [6-23]. In the present work, we restrict our study to (f3~ (3~) 2v 
decay of 94 - 96 Zr, 98 - 100 Mo, 104 Ru, 110 Pd, 128,130^ and iso Nd i sot0 pes for the 0+ 0+ transition only. 

The study of ((3~ (3~) 2v decay is quite interesting from the nuclear structure point of view. The 
((3~f3~) 2u decay has been experimentally observed in case of ten nuclei namely, 48 Ca, 76 Ge, 82 Se, 96 Zr, 
100 Mo, 116 Cd, 128 ' 130 Te, 150 Nd and 238 U out of 35 possible candidates [24] and one can extract nuclear 
transition matrix elements (NTMEs) from the observed half-lives. Using the average half-lives [28] and 
the phase space factors 141, the extracted NTMEs M 2v vary from 0.0152±0.0008 (0.0238±0.0013) to 
0.1222±0.0034 (0.1909±0.0053) for g A = 1.25 (1.00) corresponding to 130 Te and 100 Mo respectively. A 
comparison between the theoretically calculated and experimentally extracted NTMEs provides a cross- 
check on the reliability of different nuclear models used for the calculation of NTMEs. 

It is observed in all cases of {(3~ (3~) 2u decay that the NTMEs M 2u are sufficiently quenched. The 
calculation of M 2v requires the knowledge of the (3~ or (p, n) amplitude for the initial nucleus and the (3 + 
or (n, p) amplitude of the final nucleus, which in turn requires a complete set of states of the intermediate 
nucleus in addition to the initial and final nuclear states. The understanding as well as realization of 
this quenching mechanism is the main motive of all the theoretical calculations. In solving this problem, 
different nuclear models and nuclear structure scenarios have been applied. Nuclear models, which are 
used in the calculation of NTMEs of (J3~ (3~) 2v decay, can be broadly classified in to shell model and 
its variants, quasiparticle random phase approximations (QRPA) and extensions to it and alternative 
models. The details about these models -their advantages as well as limitations- have been discussed by 
Suhonen and Civitarese [TS] and Faessler and Simkovic [15] . 

The shell-model is the best choice for the calculation of the NTMEs as it attempts to solve the 
nuclear many-body problem as exactly as possible. The large scale shell model calculations by Caurier et 
al. are more realistic in which 76 Ge, 82 Se, 124 Sn, 128 Te, 130 Te and 136 Xe have been studied [25]. The M 2v 
of 82 Se is calculated exactly and those of 76 Ge and 136 Xe are dealt in a nearly exact manner [26]. In the 
QRPA model, Vogel and Zirnbauer were the first to provide an explanation of the observed suppression of 
M 2v by a proper inclusion of ground state correlations through the proton-neutron p-p interaction in the 
5=1, T=0 channel and the calculated half- lives are in close agreement with all the experimental data [27] , 
The QRPA frequently overestimates the ground state correlations and to cure the strong suppression of 
M 2v due to increase in the attractive proton-neutron interaction, several extensions of QRPA have been 
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proposed. The most important proposals are inclusion of proton-neutron pairing, renormalized QRPA, 
higher QRPA, multiple commutator method (MCM) and particle number projection. However, none of 
the above methods is free from ambiguities [19] . Alternative models, as the operator expansion method 
(OEM), the broken SU(4) symmetry, two vacua RPA, the pseudo SU(3) and the single state dominance 
hypothesis (SSDH) have their own problems [18] . 

The {(3~(3~) 0v decay has not been experimentally observed hitherto, and only limits on half- 
lives of (/3~/3~) 0l/ decay are available. Klapdor and his collaborators have reported the (fl~fl~) 0u 
decay of 76 Ge in Heidelberg-Moscow experiment [29]. However, it is felt that the reported result 
needs independent verification by other experiments [3TJ1 123] • The observed half- life limits permit to 
extract limits on various effective lepton number violating parameters, namely Majorona neutrino mass, 
coupling of left and right handed weak currents, mass of right handed heavy neutrino, mass of right 
handed W-boson, intergeneration Yukawa coupling constants of SUSY models, leptoquark coupling 
constants, compositeness scale, constraint on sterile neutrinos, Majoron coupling constants, VEP and 
VLI parameters using the theoretically calculated NTMEs of {(3~ (3~) 0v decay. In order to get accurate 
effective lepton number violating parameters, one has to calculate reliable NTMEs for ((3~ (3~) 0v decay. 
In the absence of experimental data, it is difficult to judge the reliability of wave functions involved 
in the calculation of NTMEs of (P~P~) 0l/ decay. Usually, the reliability of wave functions is tested by 
reproducing the experimentally extracted NTMEs M% v as both the modes involve same set of initial 
and final wave functions although the nuclear transition operators are sensitive to different spin-isospin 
correlations in {(3~ (3~) 2v an d (/3~f3~) Ql/ decay modes. 

The two main ingredients deciding the structure of nuclei participating in (3(3 decay are the pairing 
and deformation degrees of freedom. The crucial role of deformation on NTMEs Mi v has been predicted 
in the case of {(3~ (3~) 2v decay of 100 Mo and 150 Nd [311 E2]- The existence of an inverse correlation 
between the GT strength and quadrupole moment has been already shown by Auerbach et al [33] and 
Troltenier et al [33]. The effect of deformation on the distribution of the Gamow- Teller strength and 
/3-decay properties has been studied using a quasiparticle Tamm-Dancoff approximation (TDA) based on 
deformed Hartree-Fock (DHF) calculations with Skyrme interactions [35] and in deformed self consistent 
HF+BCS+QRPA method with Skyrme type interactions [36]. Nacher et al [37] have presented a novel 
method of deducing the deformation of N = Z nucleus 76 Sr, based on the comparison of the experimental 
GT strength distribution B(GT) from its decay with the results of QRPA calculations. A deformed 
QRPA formalism, using deformed Woods-Saxon potentials and deformed Skyrme Hartree-Fock mean 
fields, was developed to describe simultaneously the energy distributions of the single-/3 GT strength and 
the {(3~ (3~) 2u decay matrix elements [38]. The deformation effect on the {(3~ (3~) 2v decay for ground-state 
transition of 76 Ge was studied in the framework of the deformed QRPA with separable Gamow- Teller 
(GT) residual interaction [39] . 

In the light of above discussions, the Projected Hartree-Fock Bogoliubov (PHFB) model is a 
convenient choice as an alternative model in which the pairing and deformation degrees of freedom are 
incorporated on equal footing and the rotational symmetry is restored by projection technique providing 
wave functions with good angular momentum for the parent and daughter nuclei involved in (3(3 decay. 
However, the PHFB model is unable to provide information about the structure of the intermediate 
odd-odd nuclei in its present version and hence, on the single (3 decay rates and the distribution of GT 
strength. In spite of this limitation, the PHFB model, in conjunction with pairing plus quadrupolc- 
quadrupole (PQQ) interaction [40], has been successfully applied to study the {(3 ± (3 ± ) 2l) decay for the 
+ — » + transition where it was possible to describe the lowest excited states of the parent and daughter 
nuclei along with their electromagnetic transition strengths, as well as to reproduce their measured (3(3 
decay rates [UJ 35] [J3] . In the PHFB model, the role of deformation in reproducing realistic NTMEs 
M-i v has also been investigated and it has been observed that there exists an inverse correlation between 
the latter and the former [HI l42l [43] . 

In the present work, we add a hexadecapole-hexadecapole interaction term V{HH) to the standard 
PQQ interaction to check the stability of our previous results of ((3~ (3~) 2v decay with respect to the 
change in effective two-body interaction. In variation-after-projection (VAP) framework, the pairing 
plus quadrupole-quadrupole plus hexadecapole-hexadecapole (PQQHH) interaction has been successfully 
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applied to study the yrast spectra of 68 ~ 76 Ge, 72 ~ 78 Se, 74 ~ 82 Kr, 10 °-ios Zr an d ioo-io8 Mo i soto pes 
[44] . The present paper is organized as follows. In section 2, we present the theoretical formalism 
to calculate the NTME M 2v in the PHFB model in conjunction with summation method using PQQHH 
interaction. The expressions to calculate the spectroscopic properties, specifically, the yrast spectra, 
reduced B(E2:0 + — > 2+) transition probabilities, static quadrupole moments Q(2 + ) and ^-factors 
g(2 + ) are given by Dixit et al 45J. The calculated spectroscopic properties of 94,96 Zr, 94 < 96 < 98 < 100 Mo, 
98,ioo,i04 RU) io4,no pdj no Cdj 128,130^^ 128,130^^ i50 Nd and i50g m nuclei are compare d with the observed 

experimental data in section 3 and there by, we check the "goodness of wave functions" . The same wave 
functions are used to study the li~) 2v decay of 94 < 96 Zr, 98 < 100 Mo, 104 Ru, 110 Pd, i28,i30 Xe and i 5 o Nd 
nuclei for the + — > + transition and the results are presented in the same section. We also examine the 
effect of deformation on NTMEs M 2v by varying the strength of QQHH part of the effective two-body 
interaction in section 3. While presenting the theoretically calculated results using PQQHH interaction, 
we also give our previous results [41| 143] for comparison, which were calculated using PQQ interaction. 
Finally, we present some concluding remarks in section 4. 



2. Theoretical framework 



The inverse half-life of the (f3 f3 ) 2 decay for the + — » + transition [21 033 US] is given by 

[r 1 2 ; 2 (o + ^o+)]- 1 = G 2 ,|M 2 ,| 2 (i) 

where the phase space factor G 2v can be calculated with good accuracy [TQl [14] and the nuclear model 
dependent NTME M 2v is written as 

(0+|K+||l+)(l+||ar+||0+) 

N 



M 2v = J2 



N 



Eq + Eff — Ej 



with 



E 



2 ( E i _ e f) - 2*3/3/3 



(3) 



To evaluate (|5J), it is required to sum over all the intermediate 1^ states . However, it is not possible 
to study the structure of odd-odd nuclei in the present version of the PHFB model. Hence, we carry 
out the summation over the intermediate states by using the summation method given by Civitarese and 
Suhonen |46j . Using summation method, the M 2v is expressed as 



M 2v 

where Y a is given by 
and 



with 




o 



E 

A=0 



D\Y 



(4) 

(5) 
(6) 



D x Y, t = [H, [H, ,[H,Y^] ]f times ) (7) 

When the GT operator commutes with the effective two-body interaction, the ([4]) can be further simplified 
to 



M 2v = 



(0+||a-ar+T+||0+) 



(8) 
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In the case of pseudo-SU(3) model, the GT operator commutes with the two-body interaction [151 [i9l I5T)] 
and the energy denominator is a well-defined quantity without any free parameter. It has been evaluated 
exactly for {j3~ P~) 2l> HUES] an( A (e£) 2; , modes [SO] in the pseudo-SU(3) scheme. 

In the present work, we use a Hamiltonian with PQQHH type of effective two-body interaction. The 
Hamiltonian is explicitly written as 

H = H sp + V{P) + ( qq [V(QQ) + V(HH)} (9) 

where H sp denotes the single particle Hamiltonian. The pairing part of the effective two-body interaction 
V(P) is given by 

V{P) = (t) E(- 1 ) jQ+J ^ m ^ m3a « a a^a /3 (10) 

where a denotes the quantum numbers (nljm) and the state a is same as a but with the sign of m 
reversed. The QQ part of the effective interaction V{QQ) is expressed as 

V(QQ) = - (f) E £(-1) M H<Z2mI7></%2-mI*> ataj a 5 a 7 (11) 

where 



16tt 



1/2 



92„=(^J r 2 Y 2tl (9,cp) (12) 
The HH part of the effective interaction V{HH) is given as 

V{HH) =-(f) E £(-l)>|<fe|7></%4-,,|<5) a W a * a ^ ( 13 ) 

q/37(5 f 

with 

q 4l/ = r 4 r 4 „(M) (14) 

Further, C g(3 is an arbitrary parameter and the final results are obtained by setting C 99 = 1- The purpose 
of introducing it is to study the effect of deformation by varying the strength of effective two-body QQHH 
interaction. 

The model Hamiltonian used in the present work does not commute with the GT operator. Hence, 
the energy denominator is not a well-defined quantity. However, the violation of isospin symmetry for 
the QQHH part of our model Hamiltonian is negligible as will be evident from the parameters of the 
two-body interaction given in section 3. Further, the violation of isospin symmetry for the pairing part 
of the two-body interaction is presumably small. With these assumptions, the NTME M2,, of (/3~/?~) 2 „ 
decay for the + — * + transition in the PHFB model in conjunction with the summation method can 
be obtained as follows. 

In the PHFB model, states with good angular momentum J are obtained from the axially symmetric 
HFB intrinsic state |<i>o) with K—0 using the standard projection technique [FT] given by 

2 7-1-1 f 

l*oo) = J D-UmW |*o) dSl (15) 

where R(tt) and Dq (H) are the rotation operator and the rotation matrix respectively. The axially 
symmetric HFB intrinsic state |$ ) can be written as 

i$ )=n( u »»+ u ™ 6 i™^)i°) ( i6 ) 

im 

where the creation operators i>J and b\ m are defined as 

bim — E ^a,™ a <Lm EnC ^ ^im = ^^(~-*-)' + "' "^ia^da^m (17) 
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Finally, one obtains the expression for the NTME M 2i , of (/? (3 )„ decay for the + — > + transition as 

nasi] 

_ ...> j /=°i 

M 2v = J] ■ 



where 



oo 



la-ar+r+lK- ) 



„Ji=0 J/=o 
"■(Z,N) n (Z+2,N-2) 



-1/2 



n (Z,N),(Z+2,N 



(I 



ap-y8 



e>7 



(")* 



-sinOdd 



n J = 



det (l + FW/ (w)+ )] 1/a (l + i ?M / Wt )] V2 di {6)sin{6)d6 



(18) 



(19) 



and 



"•(Z,iV),(Z+2,iV-2)(0) = det (1 + Fz,Nfz+2,N-2 



1/2 



det 1+F^/^ 2)JV _ 2 (20) 



1/2 



The 7r(i/) represents the proton (neutron) of nuclei involved in the (/? /3 ) 2 decay process. The matrices 
fz,N and F z ,n{9) are given by 



fz,N — c 



i^(0) = E 



(21) 
(22) 



With the assumption that the difference in single particle energies of protons in the intermediate 
nucleus and neutrons in the parent nucleus is mainly due to the difference in Coulomb energies, one 
obtains 

A c for n v = n^,l u = U,] v = j„ 



^(^77? 1-7T1 jV) ^(j^V: jl/) 



AE s .o. splitting J " Or n v Tt^ , \ v 1>tti3i> 7^ Jtt 

where the Coulomb energy difference Ac is given by Bohr and Mottelson as [47] . 
0.70 r 



,(23) 



A 1 / 3 



(2Z + 1) 



0.76{(Z + l) 4/3 -Z 4 / 3 } 



MeV 



(24) 



The numerical calculation of M 2l/ for the nuclei involved in the (J3~ /3~) 2 „ decay involves the setup of 
matrices fz,N and Fz,n{@) given by (f2"Tj) and (|2"2"j) at 20 Gaussian quadrature points in the range (0, 7r) 
using the results of PHFB calculations, which are summarized by amplitudes (ui m , Uj m ) and expansion 
coefficients Cij <m . Subsequently, the required NTME is evaluated in a straightforward manner using (|18|) . 

It must be underlined that, in the present context, the use of the summation method goes beyond 
the closure approximation, because each proton-neutron excitation is weighted depending on its spin-flip 
or non-spin-flip character. The explicit inclusion of the spin-orbit splitting in the energy denominator, 
(|23| . implies that it cannot be factorized out of the sum in ([2]). In this sense, employing the summation 
method in conjunction with the PHFB formalism is richer than what was done in previous application 
with the pseudo SU(3) model [4"5H4"5] . 
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3. Results and discussions 

The model space and single particle energies (SPE's) are the same as in our earlier calculation on (/3~/3~ ) 2 „ 
decay for the + — > + transition [HI 03]. However, we briefly discuss in the following the model space 
and single particle energies (SPE's) used to generate the HFB wave functions for convenience. We treat 
the doubly even nucleus 76 Sr (N = Z = 38) as an inert core in case of 94 ' 96 Zr, 94,96,98,ioo MO; 98,100,104 RU) 
104,110 Pd and 110 Cd nuclei, with the valence space spanned by lpi/ 2 , 2s 1 / 2i ld 3 / 2} 1^5/2: 0g 7 / 2 , 0ff9/2 an d 
0/in/2 orbits for protons and neutrons. The lpi/2 orbit has been included in the valence space to examine 
the role of the Z = 40 proton core vis-a-vis the onset of deformation in the highly neutron rich isotopes. 
The set of single particle energies (SPE's) used here are (in MeV) e(lpi/ 2 ) = —0.8, e(0g 9 / 2 ) = 0.0, 
e(ld 5 / 2 ) = 5.4, e(2si/ 2 ) — 6.4, e(ld 3 / 2 ) — 7.9, e(0g 7 / 2 ) = 8.4 and s(0hn/ 2 ) — 8.6 for proton and neutron. 

In case of 128 < 130 Te, 128 < 130 Xe, 150 Nd and 150 Sm nuclei, we treat the doubly even nucleus 100 Sn 
(N = Z = 50) as an inert core with the valence space spanned by 2s\/ 2 , ld 3 / 2 , 1^5/2, I/7/2, O97/21 
O/ig/2 and 0/in/2 orbits for protons and neutrons. The change of model space is forced upon as in the 
model space used for mass region A w 100, the number of neutrons increase to about 40 for nuclei 
occurring in the mass region A w 130. With the increase in neutron number, the yrast energy spectra 
gets compressed due to increase in the attractive part of effective two-body interaction. The set of 
single particle energies (SPE's) used here are in MeV: e(ld 5 / 2 ) = 0.0, £(252/2) = 1-4, e(ld 3 / 2 ) = 2.0, 
£(057/2) = 4.0, e(0h 11/2 ) = 6.5 (4.8 for 150 Nd and 150 Sm), e(lf 7/2 ) = 12.0 (11.5 for 150 Nd and 150 Sm), 
£(0/19/2) = 12.5 (12.0 for 150 Nd and 150 Sm) for proton and neutron. 

The HFB wave functions are generated using an effective Hamiltonian with PQQHH type of two- 
body interaction. Explicitly, the Hamiltonian can be written as 

H = H sp + V(P) + Q qq [V{QQ) + V{HH)} (25) 

where H sp denotes the single particle Hamiltonian. The V(P), V(QQ) and V(HH) represent the pairing, 

quadrupole-quadrupole and hexadecapole-hexadecapole part of the effective two-body interaction. The 

C qq is an arbitrary parameter and the final results are obtained by setting the £ qq = 1. The purpose 

of introducing Q q is to study the role of deformation by varying the strength of QQ interaction. The 

strengths of the pairing interaction is fixed through the relation G p = 30 /A MeV and G n = 20 /A MeV, 

which are same as used by Heestand et al. [52] to explain the experimental <?(2 + ) data of some even-even 

Ge, Se, Mo, Ru, Pd, Cd and Te isotopes in Greiner's collective model [53]. For 94 Zr and 96 Zr, we have 

used G n = 18/ A and 22/^4 MeV respectively. The strengths of the pairing interaction fixed for 128 < 130 Te, 
i28,i30 Xej iso Nd and i50 gm are Gp = 35^4 MeV and Gn = 35^4 MeV ^ 

The strengths of the like particle components of the QQ interaction are taken as X2pp = X2nn = 
0.0105 MeV &~ 4 , where b is oscillator parameter. The strength of proton-neutron (pn) component of the 
QQ interaction X2pn is varied so as to obtain the spectra of considered nuclei namely 94:96 Zr, 94 > 96 > 98 > 100 ]y[o, 
98,ioo,io4 RUj io4,iio pd5 iio Cdj 128,130^ 128,130^ i 5 o Nd and i50 Sm in optimum agre ement with the 

experimental results. The theoretical spectra has been taken to be the optimum one if the excitation 
energy of the 2 + state E 2 + is reproduced as closely as possible to the experimental value. The prescribed 
set of parameters for the strength of QQ interaction are consistent with those of Arima suggested on the 
basis of an empirical analysis of effective two-body interaction [54] . 

The relative magnitudes of the parameters of the HH part of the two body interaction are calculated 
from a relation suggested by Bohr and Mottelson [SH] . According to them the approximate magnitude of 
these constants for isospin T — is given by 

^=2ATTZpi) ^A = 1,2,3,4-.. (26) 

and the parameters for the T = 1 case are approximately half of their T = counterparts. We take the 
X4 for T = 1 case as exactly half of the T = case. Using b = 1.0032A 1 / 6 , one obtains 

- v (I)0)V^ 4 

= 0.4884x2^~ 2/3 k~ 4 (27) 
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We fix X2pn through the experimentally available energy spectra for a given model space, SPE's, G p , G n 
and X2 PP - We present the values of X2 P n m table 1. All these input parameters are kept fixed to calculate 
other spectroscopic properties. Further, we have performed independent calculations for the parent and 
daughter nuclei involved in the (3(3 decay, whose deformations are in general different. 

3.1. Yrast spectra and electromagnetic properties 

The calculated values of excitation energies E 2 + , £4+ and E 6 + for all the nuclei of interest along with 
the experimental ones [55] are given in table 1. It can be seen from table that the theoretical spectra is 
more expanded in comparison to the experimental spectra for all nuclei although the agreement between 
the theoretically calculated and experimentally observed E 2 + is quite good. This can be taken care in 
conjunction with the VAP prescription [57] . However, our aim is to reproduce the properties of only low- 
lying 2 + state and hence, we do not invoke the VAP prescription, which will unnecessarily complicate 
the numerical calculation. 

The reduced transition probabilities B(E2:0 + — > 2 + ) are calculated for effective charges e e // = 
0.40, 0.50 and 0.60. The experimentally observed results [55J [59] are also given in the same table. It 
is noticed that the calculated and the observed B(E2:0 + — » 2 + ) values are in excellent agreement in 
case of 94 Zr, 94 Mo, 104 Ru and 104 Pd isotopes for e e // = 0.60. For the same e e //, the theoretically 
calculated B(E2:Q+ -> 2+) differ by 0.007 and 0.004 e 2 b 2 only from the experimental data for 100 Mo 
and 100 Ru isotopes. In case of 96 Zr, 96 Mo, 110 Cd, 128 > 130 Te, 128 > 130 Xe and 150 Nd isotopes, the calculated 
B(E2:0 + — > 2 + ) agree with experimentally observed values at e e // = 0.50. However, the calculated 
B(E2:0 + — > 2+) differ by 0.056 and 0.508 e 2 b 2 from the experimental results in case of 110 Pd and 150 Sm 
nuclei for e e// = 0.50. The calculated B{E2:0+ -> 2+) of 98 Mo and 98 Ru nuclei for e eff = 0.40 are in 
agreement with the experimental values. 

The theoretically calculated Q(2 + ), which are calculated for the same effective charges, i.e. e e // = 
0.40, 0.50 and 0.60, and the experimental Q(2+) results [60] are given in the same table 2. No 
experimental Q{2 + ) result is available for 94 ' 96 Zr and 128,130 Xe nuclei. For the same effective charge 
as used in case of B(E2:0 + — ► 2 + ), the agreement between the calculated and experimental Q(2 + ) is 
good for 104 Ru, 110 Pd and 150 Sm nuclei. The calculated Q(2+) are off by 0.10, 0.135, 0.02 and 0.021 eb 
in case of 98,100 Mo, 100 Ru and 150 Nd nuclei respectively from the experimental values. The theoretical 
Q{2 + ) results are quite off from the observed values for the rest of nuclei. 

The gyromagnetic factors g{2 + ) are calculated with gf = 1.0, g\ — 0.0 and gl = g v s = 0.60 and 
presented in table 2 along with the available experimental g(2 + ) data 60J. No experimental result is 
available for 96 Zr and 94,96 Mo. The theoretical g{2 + ) value of 94 Zr is a pathological case. The calculated 
g(2 + ) is 0.112 nm while the most recent measured value is —0.329 ± 0.015 nm [61]. The calculated and 
experimentally observed g(2 + ) results are in good agreement for 98 ' 100 Mo, 98 Ru, 104 Pd, 110 Cd and 128 Xe 
nuclei. The discrepancy between calculated and experimentally observed g{2 + ) values are 0.047, 0.015, 
0.089 and 0.028 nm only for 100 < 104 Ru, 110 Pd and 130 Xe nuclei respectively. The theoretically calculated 
and experimentally observed g{2 + ) results are off by 0.136, 0.257, 0.161 and 0.161 nm in case of 128 > 130 Te, 
150 Nd and 150 Sm nuclei respectively. 

It is observed from table 1 and table 2 that the calculated yrast spectra, reduced B(E2:0 + — » 2 + ) 
transition probabilities, static quadrupole moments Q{2 + ) and gyromagnetic factors g{2 + ) with PQQ 
and PQQHH type of effective two body interaction do not differ much. This establishes the stability of 
our previous results [41j [43] against the change in effective two-body interaction. 

3.2. Results of {(3" (3~) 2v decay 

For the calculation of {(3~ (3~) 2v decay rates, the phase space factors G 2v of 94,96 Zr, 98,100 Mo, 104 Ru, 
110 Pd, 128 > 130 Te and 150 Nd nuclei for the 0+ -^0+ transition are 2.304xl0" 21 yr" 1 , 1.927xl0~ 17 
yr~\ 9.709xl0" 29 yr~\ 9.434xl0~ 18 yr" 1 , 9.174xl0~ 21 yr~\ 3.984xl0~ 19 yr"\ 8.475xl0~ 22 yr" 1 , 
4.808xl0 -18 yr -1 and 1.189xl0~ 16 yr -1 respectively for gA= 1.25 [M]. However, it is more justified 
to use the nuclear matter value of gA around 1.0 in heavy nuclei. Hence, the experimental M 2ll as well 
as the theoretical T^Y 2 are calculated for gA = 1-0 and 1.25. In table 3, we compile all the available 
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experimental and the theoretical results along with our calculated M 2v and corresponding half-lives T^J 2 
of all the nuclei under consideration for the + — > + transition. We also present the Mi v extracted from 
the experimentally observed Tfj 2 using the given phase space factors. We present only the theoretical 
Ty 2 for those models for which no direct or indirect information about M 2u is available to us. 

The (/3 _ /3 _ ) 2 „ decay of 94 Zr for the + — » + transition has been investigated experimentally only 
by Arnold et al [52] and the reported limit is Ty 2 > 1.1 x 10 17 yr. The half- life calculated in the PHFB 
model for gA — 1-25 is 1.086 x 10 23 yr, which lies within the range given by Bobyk et al 63J . The 
calculated half-life Tfj 2 in QRPA model [65] is smaller than the presently calculated Tf" 2 by a factor of 
1.5 approximately for gA = 1.25. On the other hand, the half-life T^" 2 calculated in OEM [53] is larger 
than our PHFB model value for gA — 1-25 by a factor of 15 approximately. The predicted T^ 2 for 
g A = 1.0 in the PHFB model is 2.652 x 10 23 yr . 

In the case of 96 Zr^ 96 Mo, the theoretically calculated M 2u given by Stoica using SRPA [75] is 
too small in comparison to the NTME extracted from the experimental data. On the other hand, the 
calculated half-life T^J 2 in the OEM [64] is quite off from the observed experimentally observed result. The 
M 2v calculated by Engel et al [73] and Barabash et al [73] using QRPA are close to the experimentally 
observed lower limit of Wieser et al [57] for gA = 1.0. The Ty 2 calculated by in RQRPA using WS basis 
and AWS basis [7T] are 4.2 x 10 19 yr and 4.4 x 10 19 yr respectively and agree with the experimental T 2 ^ 2 
of Kawashima et al |69] . The Ty 2 predicted by Staudt et al [65j is in agreement with the experimental 
results of Barabash [55] and Wieser et al [57] . The experimentally observed Tfj 2 given by NEMO 3 
[66] and Wieser et al [67] are favored by the Ty 2 calculated in the PHFB model and SU(4) CTT [70] for 
gA = 1.25. The predicted half-life T 2 ^ of Bobyk et al [53] has a wide range and cover all the available 
experimental results. 

No experimental result for Tfj 2 of 98 Mo isotope is available so far. It has been studied theoretically 
in the QRPA [55], OEM [5J and SRQRPA [53]. The presently calculated Ty 2 = 6.36xl0 29 -1.552xl0 30 
yr for gA = 1-25 — 1.0 respectively is within the range given by SRQRPA model. The calculated T 2 ^ 2 in 
QRPA and OEM are larger than our predicted value for gA = 1.0 by approximately a factor of 2 and 4 
respectively. 

In the case of 100 Mo, the theoretically calculated M 2v in SRPA [75] is too small in comparison to 
the experimental M 2v . The (/3 _ /? _ ) 2 „ decay rate of 100 Mo calculated by Staudt et al in QRPA [55J 
and Hirsch et al in OEM [51] are off from the experimental T^J 2 . The calculated M 2v by Griffiths et al 
[5T] using QRPA model favors the results of INS Baksan [53] and LBL [ST] for g A = 1.0. On the other 
hand, the M 2v predicted by Engel et al [73] and Civitarese et al [55] are in agreement with the results 
of LBL [81], LBL collaboration [79], UC Irvine [78] and ITEP+INFN [77] for g A = 1.0. The NTMEs 
M 2v predicted in SU(4) CTT [70] and pseudo-SU3 using spherical occupation wave functions [3S] are nearly 
identical and close to the experimental result given by INS Baksan and ITEP+INFN for gA = 1.25. The 
same two M 2v are in agreement with the results of UC-Irvine 83J, ELEGANTS V, LBL and NEMO 3 
for g A = 1-0. Further, the NTME M 2u given in the PHFB model, QRPA [52] and pseudo-SU3 using 
deformed occupation wave functions [49] favor the results of UC-Irvine [55] . ELEGANTS V, LBL, LBL 
collaboration and ITEP+INFN for g A = 1.25. The results of SSDH [85] are in agreement with the 
experimental half-lives of LBL, LBL collaboration, UC-Irvine [75] and ITEP+INFN. The decay rate Tfj 2 
calculated by Bobyk et al [63j is in agreement with all the experimental results due to a large range of 
values. 

The {(3~P~) 2v decay of 104 Ru— ^ 104 Pd for the + — * + transition has not been experimentally 
investigated so far. The half- life Ty 2 has been calculated theoretically in QRPA [65] and OEM [64] , 
Our predicted Ty 2 for gA = 1-25 is approximately 3.7 times larger than that of calculated in QRPA and 
smaller than the half-life predicted in OEM by a factor of approximately 1.3. The predicted Tfj 2 for 
g A = 1.0 in the PHFB model is 5.73 x 10 22 yr. 

In 1952, Winter had studied the (f3~ fl~) 2l/ decay of 110 Pd isotope for the + — > + transition [57] , 
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The deduced half-life limit was Tfj 2 > 6.0 x 10 16 yr for (J3~0~) 2v decay mode and T*Y 2 > 6.0 x 10 17 
yr for all modes. Since then, no experiment has been attempted to study the f3~{3~ decay of 110 Pd. 
The (/?~/3~) 2l/ decay of 110 Pd has been studied theoretically in QRPA [55], OEM [53], SRPA [55] and 
SSDH [86j[H8]. The calculated iffo for g A = 1.25 in the present PHFB model is 1.74 x 10 20 yr, which is 
close to that of Semenov et al [88] and approximately twice of Civitarese et al [86] . On the other hand, 
the calculated half-life in SRPA is about 7 times larger than the presently calculated T 2 ^ for the same 
g A . The calculated half-life are 1.16 x 10 19 yr and 1.24 x 10 21 yr in QRPA and OEM respectively. The 
predicted half-life in the present work for g A — 1.0 is TfY 2 = 4.24 x 10 20 yr. 

In the case of 128 Te nucleus, the presently calculated M 2v is close to the NTME extracted from the 
experiments of Takaoka et al [5D], Manuel [55] and Barabash [35] for g A — 1.0 while it is close to those of 
Lin et al [5J and Manuel [5J] for g A = 1.25. The NTMEs M 2v calculated in SU(4) CTT [70], SSDH [55] and 
MCM 95j differ from the presently calculated M 2v by factor of 1.4 — 1.6 and are in agreement with the 
experimental NTMEs due to Lin et al [93] and Manuel [94] for g A = 1.0. The NTME M 2v calculated in 
SSDH [55] is close to the experimental M 2v of Takaoka et al [90 , Bernatovicz et al [5T] and Barabash [55J 
for g A — 1.25. The presently calculated M 2v is smaller than those calculated in QRPA [73] and WCSM 
[5J by factor of approximately 2.2 and 3.6 respectively while larger by a factor of approximately 5.5 and 
2.5 than the NTMEs calculated in SRPA [S5J and SSDH [55] respectively. The half-life lfj 2 calculated 
in QRPA [65] is close to the experimental Tfj 2 due to Lin et al [93] and Barabash [28]. The calculated 
half-lives T 2 ^ by Hirsch et al in OEM [53] and Scholten et al in IBM [55] are quite small. In the PHFB 
model, the predicted half-life Ty 2 is (1.10 - 2.68) x 10 24 yr for g A = (1.25 - 1.0). 

In the PHFB model, the calculated NTME M 2v for the (J}~P~) 2u decay of 130 Te is in agreement with 
the NTME extracted from the Milano+INFN experiment [57] for g A = 1.0, theoretical M 2v calculated 
in the shell model [25] and MCM [95]. The calculated NTME M 2v in SRPA [89] is in agreement with 
the experimental M 2u J93][90J[97] for g A = 1.25. The presently calculated M 2l/ is smaller by a factor of 
approximately 1.5 than those calculated in SU(4) CTT [70] and QRPA [74 . The M 2l/ calculated in RQRPA 
[71] is close to the experimental NTME of Bernatovicz et al [91] for g A = 1.25. The calculated half- lives 
Tfy in OEM [53], IBM [55J and WCSM 9 are quite small while the half-life Tfy calculated in SU(4) (TT 
[98] is close to the experimentally observed half-lives [90J [93] [97] . The predicted half-life Ty 2 in the 
PHFB model is (2.10 - 5.13) x 10 20 yr for g A = (1.25 - 1.0). 

In the case of 150 Nd isotope, the calculated M 2v in the PHFB model is in agreement with the 
extracted NTME from NEMO 3 [55] project for g A = 1.25. The NTME M 2v calculated in pseudo-SU(3) 
[48] is close to the experimental M 2v of NEMO 3 [55] for g A = 1.0. The calculated M 2v in SU(4) CTT 
[TP] is larger by a factor of 2.4 than the present value. The half-life T 2 ^ calculated in OEM [53] is 
close to the experimental result of ITEP+INR project while the calculated Ty 2 in QRPA 65] favors 
the UCI experiment [78]. The calculated half-life jfY 2 in the PHFB model is (1.19 - 2.92) x 10 19 yr for 
g A = (1.25 -1.0). 

In all the cases discussed in table 3, it turns out that the NTMEs M 2v calculated with PQQHH 
and PQQ effective two-body effective interactions are very close. Further, it is noticed that although 
the calculated spectroscopic properties in the PHFB model for nuclei in the vibrational limit are not in 
agreement with the experimental data as expected, the calculated NTME M 2v of 96 Zr for example is in 
good agreement with the experimentally observed data. 

3.3. Deformation effect 

To understand the role of deformation on the NTMEs M 2l/ , we investigate its variation with respect to 
the change in strength of the QQHH interaction ( qq . The results are presented in figure 1. In all cases, 
it is observed that NTMEs M 2v remain almost constant as the strength of ( qq is changed from 0.0 to 
0.6 except in case of 94 Zr and 100 Mo isotopes, where M 2u remains almost constant for the variation of 
C, qq from 0.0 to 0.8. As C, qq is further increased up to 1.5, the NTME M 2v starts decreasing except a few 
anomalies. The experimental M 2v is available for {fi~ p~) 2v decay of 96 Zr, 100 Mo, 128 > 130 Te and 150 Nd 
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isotopes. It is noteworthy that the M 2v tends to be realistic as £ qq acquires a physical value around 1.0. 
This suggest that the deformations of the HFB intrinsic states play an important role in reproducing a 
realistic M 2v . Further, it is observed in general that there is an anti-correlation between the NTME M 2v 
and the deformation parameter 13 2 . 




Figure 1. Dependence of M2„ on the strength of the QQHH interaction C, qq . 



To quantify the deformation effect on M 2 „, we define a quantity D 2v as the ratio of M 2v at zero 
deformation (( qq = 0) and full deformation (( qq = 1). This ratio of deformation effect D 2v is given by 

n _ M 2lJ (( qq = 0) 

° 2v M 2 „(C W = 1) (28) 

The values of D 2v for 94 < 96 Zr, 98 ' 100 Mo, 104 Ru, 110 Pd, 128 < 130 Te and 150 Nd nuclei for PQQHH(PQQ) type 
of effective two-body interaction are 2.23(2.29), 3.38(3.70), 1.75(1.86), 2.09(2.33), 5.19(5.47), 3.30(3.14), 
3.97(4.26), 3.31(2.89) and 6.46(5.94) respectively. These values of D 2v suggest that the M 2v is quenched 
by a factor of approximately 2 to 6.5 in the mass range A = 90— 150 due to deformation effects. Further, 
the quenching is of almost same magnitude in the calculation using both type of effective two-body 
interaction. We like to mention here that in view of the limitations of the PHFB model as mentioned, 
the deformation effect for nuclei like 94 Zr and 96 Zr may be taken as a conservative estmate. 
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4. Conclusions 

To summarize, we employ the PHFB model using PQQHH type of effective two-body interaction to 
construct the yrast band wave functions of (3(3 emitters. The overall agreement between the calculated 
and observed yrast spectra as well as electromagnetic properties of the nuclei suggests that the PHFB 
wave functions generated by fixing X2pn to reproduce the E 2 + are quite reliable. Subsequently, we employ 
the same wave functions to study the (/T /3~) 2u decay of 94 - 96 Zr, 98 > 100 Mo, 104 Ru, 110 Pd, 12 ^ 130 Te and 
150 Nd nuclei for the + — > + transition. The theoretically calculated NTMEs M 2v are compared with 
those extracted from experimentally observed half-lives Ty 2 . The agreement between the theoretically 
calculated results and experimentally observed data is quite satisfactory. We also examine the effect of 
deformation on NTMEs M 2v by varying the strength of QQHH part of the effective two-body interaction. 
It is noticed that the desired quenching of M 2v is achieved in the PHFB model through the subtle interplay 
of pairing and deformation degrees of freedom. Further, the results calculated with PQQ and PQQHH 
type of effective two-body interactions are quite similar, which underlines the stability of our previous 
results [13113]. 

Further, a reasonable agreement between the calculated and observed spectroscopic properties and 
A/2!/ of (/3~/3~) 2j , decay for the + — > + transition of considered nuclei makes us confident to employ 
the same PHFB wave functions to study the (/3~/3~) „ decay of the same nuclei, which will be reported 
soon. 

This work was partially supported by DAE-BRNS, India vide sanction No. 2003/37/14/BRNS/669, 
by Conacyt-Mexico and DGAPA-UNAM. 
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Table 1: Excitation energies Ej* (MeV) of J*" = 2+, 4+ and 6 + yrast states of 

94,96^ 94,96,98,100 MO; 98,100,104^ 104,110p d; 110^ 128,130^ 128430 X ^ 150 Nd and 150 Sm M 

Nucleus PQQ PQQHH Exp.[56, Nucleus PQQ PQQHH Exp. 56] 
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E 6 + 


9 


.7261 


9. 


.3686 








E 6 + 


3 


,5776 


3. 


,5300 


2 


.44064 


98 Mo 


X2pn 


0. 


.01955 


0. 


.01955 






98 Ru 


X2pn 





.02763 


0. 


,02649 








E 2 + 


0. 


.7892 


0. 


.7887 


0. 


.78742 




E 2 + 





.6513 


0. 


,6522 





.65241 




77 

-C/4+ 




n c oo 
.9522 


1 


nn 1 a 
.9916 


1 


.51013 




TP 


1 


n a on 
943U 


1. 


n a 7n 
,94/0 


1 


.39 (0 




E 6 + 


3. 


.3098 


3. 


.4051 


2 


.3438 






3 


.6548 


3. 


,6703 


2 


.2227 


100 Mo 


X2pn 


0. 


.01906 


0. 


.01876 






100 Ru 


X2pn 





.01838 


0. 


,01831 








E 2 + 


0. 


.5356 


0. 


.5357 


0. 


.5355 




E 2 + 





.5395 


0. 


,5402 





.5396 




TP 

-C/4+ 


1. 


a 71 n 
.4/ 19 


1. 


.4/66 


A 

1. 


1 o c n 
.1359 




TP 

-C/4+ 


1 

1 


.5591 


1 

1. 


,5847 


1 

1 


00c c 
.2265 




E e + 


2. 


.6738 


2. 


.6893 








E 6 + 


2 


.8940 


2. 


,9629 


2 


.0777 


104 Ru 


X2pn 


0. 


.02110 


0. 


.02053 






104p d 


X2pn 





.01486 


0. 


,01507 








E 2 + 


0. 


.3580 


0. 


.3578 





.35799 




E 2 + 





,5552 


0. 


,5560 





,55579 




77 

-C/4+ 


1 


1 on 
.1339 


1 


.1385 


n 

u. 


.8885 




TP 


1 


c Ton 
.5729 


1 

1. 


,6138 


1 


o oo c n 
.32359 






2 


.2280 


2 


.2486 


1 


.5563 






2 


.8790 


2. 


,9954 


2 


.2498 


110p d 


X2pn 


0. 


.01417 


0. 


.01393 






110 Cd 


X2pn 





.01412 


0. 


,01414 








E 2 + 


0. 


.3737 


0. 


.3738 


0. 


.3738 




E 2 + 





.6576 


0. 


,6585 





.6577 




TP 

£/4 + 


1. 


.1563 


1. 


1 c 
.1583 


n 

u. 


nono 
.9208 




TP 

-C/4+ 


1 

1 


7nn 
.S/U9 


1 

1. 


ono 1 
,8921 


1 

1 


.5424 




E 6 + 


2. 


.2254 


2. 


.2359 


1, 


.5739 




E 6 + 


3 


.3865 


3. 


,4728 


2 


.4799 


128 Te 


X2pn 


0. 


.02715 


0. 


.02692 






128 Xe 


X2pn 





.0360 


0. 


,02662 








E 2+ 


0. 


.7436 


0. 


.7435 


0. 


.7432 




E 2+ 





.4511 


0. 


,4420 





.4429 




E4+ 


2 


.0458 


2 


.0130 


1 


.4971 




-E4+ 


1 


.4263 


1. 


,3444 


1 


,0329 




E 6 + 


3. 


.7363 


3. 


.6647 


1 


.8111 




E e + 


2 


.7976 


2. 


,5404 


1 


,7370 


130 Te 


X2pn 


0. 


.01801 


0. 


.01890 






130 Xe 


X2pn 





.02454 


0. 


,02281 








E 2 + 


0. 


.8393 


0. 


.8395 


0. 


.8395 




E 2+ 





.5385 


0. 


,5384 





,5361 




-E4+ 


1, 


.7741 


1, 


.8085 


1. 


.6325 




E4+ 


1 


.5496 


1. 


,5268 


1 


,2046 




-£'6+ 


3. 


.0833 


3. 


.1283 


1. 


.8145 




E 6 + 


2 


.7831 


2. 


,7400 


1 


,9444 


150 Nd 


X2pn 


0. 


.02160 


0. 


.02228 






150 Sm 


X2pn 





.01745 


0. 


,01730 








E 2+ 


0. 


.1307 


0. 


.1300 


0. 


.13012 




E 2+ 





.3328 


0. 


,3359 





,33395 




E 4 + 


0. 


.4320 


0. 


.4305 


0. 


.3815 




Eaa- 


1 


.0156 


1. 


,0290 





,77335 




E 6 + 


0. 


.8960 


0. 


.8958 








Eq + 


1 


.9185 


1. 


,9504 


1 


,27885 
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Table 2: Comparison of calculated and experimentally observed reduced transition probabilities 
B(E2:0+ -> 2+), static quadrupole moments Q(2+) and g-factors g{2+) of 94 ^ 96 Zr, 94,96,98,ioo MOj 

98,100,104 RU) 104,H0p d) 110^ 128,130^ 128,130 X ^ 150 Nd and 150 Sm nudeL Here and Q (2+) 

are calculated for effective charge e p = 1 + e e / / and e n = e e t / . Here f and \ denote the calculation with 
PQQHH and PQQ type of effective two-body interaction respectively. 
Nucleus B{E2:0 + -> 2+) (e z b' z ) Q(2+) (eb) g(2+)~^mT) 



Theory Exp. [58] Theory Exp. [60] Theory Exp. [60] 

e e // e e ff 

0.40 0.50 0.60 0.40 0.50 0.60 



94 



Zr 



0.056 
0.046 



0.075 
0.062 



0.097 
0.081 



0.066±0.014* 

0.081±0.017 

0.056±0.014 



-0.207 
-0.168 



-0.240 
-0.195 



-0.272 
-0.222 



f 0.112 
$ 0.121 



-0.329±0.015 Q 

-0.26±0.06 

-0.05±0.05 



04 



Mo 



0.148 
0.148 



0.187 
0.188 



0.232 
0.232 



0.203±0.004* 

0.230±0.040 

0.270±0.035 



f -0.347 -0.391 -0.435 -0.13±0.08 f °-343 
t -0.347 -0.391 -0.435 j 0.343 



(Ki 



Zr 



!)<) 



Mo 



f 0.046 0.063 
j 0.044 0.060 



0.250 
0.265 



0.317 
0.335 



0.081 
0.078 

0.391 
0.413 



0.055±0.022* f 

i 



0.271±0.005* 

0.310±0.047 

0.302±0.039 



-0.009 
-0.012 

-0.453 
-0.466 



-0.011 
-0.015 

-0.509 
-0.524 



-0.014 
-0.018 

-0.566 
-0.582 



-0.20±0j 



t 0.297 

t 0.254 

t 0.535 

t 0.563 



0.246 
0.234 



0.316 
0.302 



0.395 
0.378 



0.267±0.009* 

0.259±0.010 

0.260±0.040 



f -0.450 -0.510 -0.570 -0.26±0.09 f °- 405 
j -0.439 -0.498 -0.557 t 0.376 



0.34±0.18 



^Ru 



0.389 
0.433 



0.488 
0.543 



0.599 
0.665 



0.392±0.012* 

0.411±0.035 

0.475±0.038 



f -0.565 -0.633 -0.700 -0.20±0.09 f 0.510 
t -0.596 -0.667 -0.739 -0.03±0.14 j 0.528 



0.39±0.30 



100 



Mo 



0.306 
0.320 



0.394 
0.412 



0.493 
0.515 



0.516±0.010* 

0.511±0.009 

0.526±0.026 



f -0.500 -0.567 -0.635 -0.42±0.09 f 0.467 
$ -0.512 -0.581 -0.650 -0.39±0.08 j 0.477 



0.34±0.18 



100 



Ru 



0.306 
0.308 



0.390 
0.393 



0.484 
0.488 



0.490±0.005* 

0.493±0.003 

0.494±0.006 



-0.501 
-0.503 



-0.565 
-0.568 



-0.630 
-0.633 



-0.54±0.07 
-0.40±0.12 
-0.43±0.07 



0.363 
0.355 



0.47±0.06 
0.51±0.07 



104 



Ru 



0.558 
0.572 



0.714 
0.732 



0.890 
0.912 



0.820±0.012* 

0.93±0.06 

1.04±0.16 



-0.676 
-0.684 



-0.765 
-0.774 



-0.854 
-0.864 



-0.76±0.19 
-0.70±0.08 
-0.66±0.05 



0.345 
0.339 



0.41±0.05 
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Nucleus B{E2:Q+ -<• 2+) {e z h z ) Q(2+) (eb) g(2+) (nm) 



Theory Exp. [58] Theory Exp. [60] Theory Exp. [60] 
e e // e eff 

0.40 0.50 0.60 0.40 0.50 0.60 



104 



Pd 



0.371 
0.361 



0.473 
0.460 



0.586 
0.571 



0.535±0.035* 

0.61±0.09 

0.535±0.035 



f -0.550 -0.621 -0.692 -0.47±0.10 f 0.458 
t -0.543 -0.613 -0.682 $ 0.439 



0.46±0.04 
0.40±0.05 
0.38±0.04 



no 



Pd 



0.471 
0.479 



0.604 
0.614 



0.754 
0.766 



0.870±0.040* 

0.780±0.120 

0.820±0.080 



-0.620 
-0.626 



-0.702 
-0.708 



-0.784 
-0.791 



-0.72±0.14 f 
-0.55±0.08 $ 
-0.47±0.03 



0.489 
0.478 



0.37±0.03 
0.35±0.03 



no 



Cd f 0.407 0.522 0.653 
| 0.427 0.548 0.685 



0.450±0.020* 

0.504±0.040 

0.467±0.019 



-0.575 
-0.590 



-0.651 
-0.668 



-0.728 
-0.746 



-0.40±0.04 f 
-0.39±0.05 t 
-0.36±0.08 



0.377 
0.358 



0.31±0.07 
0.285±0.055 



12s 



IV f 0.300 0.384 0.479 
% 0.298 0.381 0.474 



0.383±0.006* 

0.380±0.009 

0.378±0.007 



f -0.498 -0.563 -0.628 -0.14±0.12 f 0.526 0.35±0.04 
X -0.496 -0.561 -0.626 -0.06±0.05 % 0.514 0.31±0.04 



12* 



Xe f 0.571 0.729 0.906 
| 0.637 0.819 1.024 



0.750±0.040* 

0.790±0.040 

0.890±0.230 



-0.685 
-0.724 



-0.774 
-0.820 



-0.862 
-0.917 



f 0.439 
X 0.400 



0.41±0.07 
0.31±0.03 



130 



IV f 0.242 0.304 0.373 
X 0.231 0.289 0.354 



0.295±0.007* 

0.290±0.011 

0.260±0.050 



f -0.448 -0.502 -0.556 -0.15±0.10 f 0.667 0.33±0.08 
X -0.438 -0.490 -0.542 % 0.679 0.29±0.06 



130 



Xe f 0.478 0.60 5 0.746 
| 0.493 0.624 0.769 



0.65±0.05* 

0.631±0.048 

0.640±0.160 



-0.627 
-0.637 



-0.705 
-0.716 



-0.783 
-0.795 



f 0.478 
X 0.463 



0.38±0.07 
0.31±0.04 



150 



Nd 



f 2.172 
| 2.132 



2.632 3.136 
2.580 3.070 



2.760±0.040* 

2.640±0.080 

2.670±0.100 



f -1.335 -1.469 -1.604 -2.00±0.51 f 0.622 0.422±0.039 
X -1.322 -1.455 -1.587 % 0.636 0.322±0.009 



150 



Sm f 1.741 2.098 2.488 
X 1.707 2.056 2.437 



1.350±0.030* 

1.470±0.090 

1.440±0.150 



f -1.193 -1.310 -1.426 -1.32±0.19 f 0.604 0.385±0.027 
X -1.182 -1.297 -1.412 -1.25±0.20 \ 0.592 0.411±0.032 



* Average B(E2) values 59 , a [SI] 
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Table 3: Experimental half-life T 2 ^ and corresponding NTME M 2v for the 0+ 



98,100 



Mo, 



'Ru, 



3 Pd. 



1/2 

128,130 Te and 



0+ transition of 94 - 96 Zr, 
150 Nd nuclei along with theoretically calculated NTME M 2 „ and 



half-life T?f n in different nuclear models. The numbers corresponding to (a) and (b) are calculated for 



L l/2 

gA — 1-25 and 1.0 respectively. 
calculation. 



"gch." denotes the geochemical experiment and "*" denotes the present 



Model 


Ref. 


Theory 
\M 2v \ 




rp2l/ 

J l/2 


PHFB* 


PQQHH 


0.063 


(a) 


10.86 








(b) 


26.52 




PQQ 


0.076 


(a) 


7.51 








(b) 


18.34 


SRQRPA 


m 






3.08-659 


OEM 


M 






168 


QRPA 


m 






6.93 


PHFB* 


PQQHH 


0.059 


(a) 


1.47 








(b) 


3.59 




PQQ 


0.058 


(a) 


1.56 








(b) 


3.80 












SU(4) ffT 


m 


0.0678 


(a) 


1.13 








(b) 


2.76 


RQRPA 1, 


eh 






4.4 


RQRPA* 


ED 






4.2 


QRPA 


EH 


0.12-0.31 


(a) 


0.054-0.36 








(b) 


0.13-0.88 


SRPA 


m 


0.022 


(a) 


10.72 








(b) 


26.18 


OEM 


M 






20.2 


QRPA 


m 






1.08 


QRPA 


ED 


0.124 


(a) 


0.34 








(b) 


0.82 


PHFB* 


PQQHH 


0.127 


(a) 


6.36 








(b) 


15.52 




PQQ 


0.130 


(a) 


6.09 








(b) 


14.87 


SRQRPA 


M 






4.06-15.2 


OEM 


m 






61.6 


QRPA 








29.6 


PHFB* 


PQQHH 


0.104 


(a) 


9.71 








(b) 


23.72 




PQQ 


0.104 


(a) 


9.79 








(b) 


23.90 


SSDH 


m 




(a) 


7.15-8.97 


SRQRPA 








5.04-16800 



Nuclei 



Project 



Experiment 



\M; 



2v\ 



94 Zr NEMO 2 

(10 22 yr) 



96 Zr 

(10 19 yr) 



98 Mo 
(10 29 yr) 



62 



>l.lxl0" 



(a) <62.815 

(b) <98.148 



NEMO 3 


[66] 


2.0±0.3±0.2 


(a) 


0.051. 








(b) 


0.080 


gch. 


[67] 


0.94±0.32 


(a) 


0.074; 








(b) 


0.116. 


NEMO 2 


m 


2.0^j±0.5 


(a.) 


o.o5i; 








(b) 


0.080; 


gch. 


m 


3.9±0.9 


(a) 


0.036. 








(b) 


0.057; 


Average 


m 


, 4 +3.5 


(a) 


0.061. 


Value 






(b) 


0.095; 


Recommended 


[28] 


2.0±0.3 


(a) 


o.o5i; 


Value 






(b) 


0.078; 



100 Mo NEMO 3 
(10 18 yr) 

NEMO 3 
gch. 



Eg 
m 



7.11±0.02±0.54 



7.72±0.02±0.54 



2.1±0.3 



(a) 

(b) 
(a.) 
(b) 
(a) 
(b) 



+0.0079 
-0.0054 
+0.0123 
-0.0084 
+0.0172 
-0.0101 
+0.0269 
-0.0158 
+0.021 
-0.012 
+0.033 
-0.019 
+0.0051 
-0.0036 
+0.0080 
-0.0056 
+0.0151 
-0.0283 
+0.0235 
-0.0443 
+0.0043 
-0.0034 
+0.0067 
-0.0054 



n 1oo +0.0051 
U.1ZZ_ 0045 

n , q , +0.0080 
u - lyl -0.0071 
n ii 7+0.0045 
U- 11 ' -0.0040 

n io^+O- 0070 
U.lOJ_g 0063 
n 99 r+0.0180 
U.zzo_ 0145 
n ori +0.0281 
U.O01_ n 0227 



Influence of the hexadecapole deformation on the two neutrino double-^ decay 
Table 3 continued 



19 



Nuclei 



Project 



Experiment 



Theory 



\M, 



Model 



Ref. 



I Mo 



rp2l/ 

1 l/2 



ITEP+INFN 

UC Irvine 

LBL+MHC+ 
UNM+INEL 
NEMO 



mi 
m 
m 
m 



LBL |ST] 
ELEGANTS V [82] 



UC Irvine 

INS Baksan 

Average 
Value 
Average 
Value 



83. 

ED 

m 



104 Ru 
(10 22 yr) 



no pd 
(10 20 yr) 



87 



7.2±0.9±1.8 (a) 

(b) 



6.82 



+ g;l±0.68 (a) 



7-et 2 : 2 



9.5±0.4±0.9 
9.7±4.9 



11.5 



+3.0 
-2.0 



11.6 



+3.4 



3.3 



+2.0 
-1.0 



8.0±0.6 
7.1±0.4 



>6.0xl0" 



(b) 
(a) 
(b) 
(a) 
(b) 
(a) 
(b) 
(a) 

(b) 
(a) 

(b) 
(a) 

(b) 
(a) 
(b) 
(a) 
(b) 



O.mijj; 
0.190l°' 
0.125±g 
0.1951°, 

o.nsio 

0.1851°, 
0.106l°' 
0.1651°, 
0.1051°, 
0.1631°,; 
0.096l°' 
0.150lo 



0.0961°,' 
0.149±n 
0.1791° 
0.2801° 

0.1151°; 

0.1801°' 
0.1221° 

0.1911°; 



0321 
0179 

0279 
0128 
0087 
0200 
0136 
0127 
0141 
0198 
0220 
008 

M"7 

013 

010 

0441 

0193 

0689 

0302 

0096 

0105 

0150 

0164 

0035 

0115 

0051 

0180 

0355 

0378 

0551 

0.591 

0054 

0047 

0081 

0074 

0036 

0033 

0050 

0052 



(a) <6.468 

(b) < 10. 106 



SU(4) CTT 


M 


0.1606 


(a) 


4.11 








(b) 


10.03 


SSDH 


m 


0.18 


(a) 


3.27 








(b) 


7.99 


SRPA 


m 


0.059 


(a) 


30.45 








(b) 


74.34 


P SU(3)+ 


m 


0.152 


(a) 


4.59 








fb) 


11.2 


P SU(3)++ 


m 


0.108 


(a) 


9.09 








(b) 


22.19 


OEM 


m 






35.8 


QRPA 


m 


0.101 


(a) 


10.39 








(b) 


25.37 


QRPA 


[SB 


0.256 


(a) 


1.62 








(b) 


3.95 


QRPA 








1.13 


QRPA 




0.211 


(a) 


2.38 








(b) 


5.81 


PHFB* 


PQQHH 


0.068 


(a) 


2.35 








(b) 


5.73 




PQQ 


0.068 


(a) 


2.35 








(b) 


5.73 


OEM 


164] 






3.09 


QRPA 


m 






0.629 


PHFB* 


PQQHH 


0.120 


(a) 


1.74 








(b) 


4.24 




PQQ 


0.133 


(a) 


1.41 








(b) 


3.44 


SSDH 


\m 






1.6 


SSDH 


[86] 


0.19 


(a) 


0.7 








(b) 


1.70 


SRPA 


[89] 


0.046 


(a) 


11.86 








(b) 


28.96 


OEM 


M 






12.4 


QRPA 


m 






0.116 
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Nuclei 



Project 



Experiment 



Theory 



IM, 



Model 


Kel. 






J l/2 


PHFB* 


PQQHH 


0.033 


(a) 


1.10 








0) 


2.68 




PQQ 


0.033 


(a) 


1.05 








(b) 


2.55 


SSDH 


88 


0.048 


(a) 


0.51 








(b) 


1.29 


SM 








0.5 


SU(4) CTT 


M 


0.053 


(a) 


0.42 








(b) 


1.06 


SSDH 


[86] 


0.013 


(a) 


6.98 








(b) 


17.65 


MCM 


m 


0.046 


(a) 


0.56 








fbl 


1.41 


SRPA 




0.006 


(a) 


32.78 








(b) 


82.87 


OEM 


ED 






0.21 


QRPA 


EH 






2.63 


QRPA 


HH 


0.074 


(a) 


0.22 








(b) 


0.54 


IBM 


[96] 






0.09 


WCSM 


I7T1 

.9 


0.120 


(a) 


0.08 








(b) 


0.21 


PHFB* 


PQQHH 


0.031 


(a) 


2.10 








(b) 


5.13 




PQQ 


0.042 


(a) 


1.16 








(b) 


2.82 


SM 


M 


0.030 


(a) 


2.3 








(b) 


5.84 


SU(4) aT 


m 


0.0468 


(a) 


0.95 








(b) 


2.32 


RQRPA 1, 


m 


0.009 


(a) 


25.68 








(b) 


62.70 


RQRPA* 


m 


0.009 


(a) 


25.68 








(b) 


62.70 


MCM 


m 


0.028 


(a) 


2.65 








(b) 


6.48 


SU(4) CTT 


[98] 






7.0 


SRPA 


[89] 


0.016 


(a) 


8.12 








(b) 


19.84 


OEM 


[64] 






0.79 


QRPA 


[65] 






18.4 


QRPA 


m 


0.049 


(a) 


0.87 








(b) 


2.12 


IBM 


m 






0.17 


WCSM 


m 


0.114 


(a) 


0.16 








(b) 


0.40 



12 s 



Te 
(10 24 yr) 



130 Te 

(10 20 



yr) 



gch. 
gch. 
gch. 
gch. 
gch. 



90] (2.2±0.3) 

H] 7.7±0.4 

92] 2.0 

93] 1.8±0.7 

94] 1.4±0.4 

22] 7.2±0.3 



Average 
Value 

Recommended [28] 2.5±0.4 
Value 



Milano+INFN 

gch. 

gch. 

gch 

Average 
.Value 

Recommended 
Value 



97 



90 



28 



6.1±1.4t 2 ; 9 



7.9±1.0 



En 27.o±i.o 



7.5±0.3 



1221 27±1.0 



9.0±1.0 



(a) 

(b) 

(a) 

(b) 

(a) 

(b) 

(a) 

(b) 

(a) 

(b) 

(a) 

(b) 

(a) 

(b) 



(a) 
(b) 

(a) 
(b) 

(a) 
(b) 
(a) 
(b) 

(a) 
(b) 
(a) 
(b) 



n0 o+. 00018 
uz,3 -0.0014 
n o R +0.0027 
u,3D -0.0022 
mo +0.0003 
ulz -0.0003 
+0.0005 



.019 
.024 
.038 
.026±o 

,+0. 



-0.0005 



.040; 

.029±o 

.0451°; 
,+o. 



.0131°; 
.0201°,; 
.0221°; 
.0341°; 



0071 
0039 
0111 
0060 
00.-,:', 
0031 
l H 'S3 
0053 

( 3, 

0003 
0004 
0004 
0014 
0012 
0022 
0019 



n ma+ 0232 
U.UIS^q 0043 

n29+° 0362 
u.uzy_ Q 006g 

n mc+0 0011 
°- 016 -0.0009 



0.025] 
0.009^ 
0.014^ 



f0.0018 
-0.0015 
f0.0002 
-0.0002 
,,+0.0003 
-0.0002 
n n1 7 +0.0003 
u - ul ' -0.0003 
026+ - 0005 

U.UZD_ Q q 005 

n09 + °- 0002 
u.uuy_ 0002 

nnl4 +o.ooo3 

U.U14_Q QQQ2 

015 +00 °° 9 
u.uio_ 0008 

024+ 0014 
U - UZ4 -0.0012 
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Table 3 continued 



Nuclei 



Project 



Ref. 



Experiment 

1 l/2 



Theory 



I Ma 



Model 


r? ^.r 
Ket. 


1 A/T 1 




T l/2 


PHFB* 


PQQHH 


0.027 


(a) 


11.95 








(b) 


29.17 




PQQ 


0.033 


(a) 


7.89 








(b) 


19.27 


SU(4) CTT 




0.0642 


(a) 


2.04 








(b) 


4.98 


pSU(3) 


\m 


0.055 


(a) 


2.78 








(b) 


6.79 


OEM 


\m 






16.6 


QRPA 


m 






7.37 



150 Nd NEMO 3 
(10 18 yr) 

UCI 

ITEP +INR 

ITEP +INR 

Average 
Value 
Average 
Value 



66 



99 



22 



9.7±0.7±1.0 



[78] G.75to.ll ± °- 68 



18. 



3±|;| ± 1.9 



100 17 



7.0 



hio 

-5.0 



±3.5 



+ 11.8 
-0.3 



7.8±0.7 



(a.) 
(b) 



(a) 
(b) 
(a) 
(b) 
(a) 

(b) 
(a) 

(b) 
(a) 

(b) 



029+ - 0030 
u.uza_ 0023 

n n/ic+0-0047 
°- 046 -0.0036 



U.IMO_ Q25 

055+ 00051 
u.uoo_ 0038 

021 +0.0043 
u.uzi_ 0036 

n n oo+0.0067 
u.uoo_ 0056 

n n9o +0.0092 
u.uzz_ 0056 

n35+ - 0144 

n n ot;+0.0008 

U.UOO_Q Q 135 

054+0.0012 

U.UD4_ Q211 
n nqQ +0.0016 

u.u.m_ 0014 
051+ 0025 

u - uol -0.0022 



^AWS basis; ^WS basis; + Sphcrical occupation wave functions; ++ Dcformed occupation wave functions 



